Enhanced coherent dynamics near a transition between neutral quantum-paraelectric 
and ionic ferroelectric phases in the quantum Blume-Emery-Griffiths model 



Kenji Yonemitsu 1, 2, 



00 

o 
o 

(N 
•*-> 
o 

O 



1 Institute for Molecular Science, Okazaki 444-8585, Japan 
2 Department of Functional Molecular Science, Graduate University for Advanced Studies, Okazaki 444-8585, Japan 

(Dated: October 11, 2008) 

Nonequilibrium dynamics are studied near the quantum phase transition point in the one- 
dimensional quantum Blume-Emery-Griffiths model. Its pseudo-spin component S z represents an 
electric polarization, and (S z ) 2 corresponds to ionicity, in mixed-stack charge-transfer complexes 
that exhibit a transition between neutral quantum-paraelectric and ionic ferroelectric (or antifer- 
roelectric) phases. The time-dependent Schrodinger equation is solved for the exact many-body 
wave function in the quantum paraelectric phase. After impact force is introduced on a polarization 
locally in space and time, polarizations and ionicity coherently oscillate. The oscillation amplitudes 
are large near the quantum phase transition point. The energy supplied by the impact flows linearly 
into these oscillations, so that the nonequilibrium behavior is uncooperative. 
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I. INTRODUCTION 

Correlated electron systems show various intriguing 
properties, which are not only scientifically interesting 
but potentially important for future applications. Their 
nonequilibrium dynamics have much room for searching 
novel functions. Among them, photoinduced phase tran- 
sitions attract much attention because macroscopic prop- 
erties are drastically changed by a cooperative effect often 
on an ultrafast time scale and in a controlled manner £££ 

During photoinduced dynamics in quasi-one- 
dimensional electron systems, coherent oscillations are 
often observed, e.g., in the mixed-stack charge-transfer 
tetrathiafulvalene-p-chloranil (TTF-CA) complex (in 
both short-time^ and long-time^ behaviors) , the quarter- 
filled-band charge-ordered organic (EDO-TTF) 2 PF 6 
(EDO-TTF=ethylenedioxy-tetrathiafulvalene) salt^ 
and potassium tetracyanoquinodimethane (K-TCNQ) 
crystals^ and their origins are theoretically investigated 
in itinerant electron models for TTF-CA (with short- 
time^ and long-time^ behaviors), (EDO-TTF) 2 PF 6 , 10 
and K-TCNQ;— In most of the cases, electron-phonon 
interactions are substantial. In theoretical works, 
phonons are often treated as classical variables in order 
to reduce the numerical burden when the electronic 
degrees of freedom are taken into account. However, 
quantum nature is indispensable to maintaining coherent 
oscillations at not-too-low temperatures where the clas- 
sical Maxwell distribution of thermal lattice fluctuations 
would considerably deteriorate the coherence. 

Here, we focus on the quantum phase transition re- 
alized in mixed-stack charge-transfer complexes com- 
posed of 4,4'-dimethyltetrathiafulvalene (DMTTF) and 
tetrahalo-p-benzoquinones (QBr n Ci4-„) i 12 i 13 For n—2, 
3, and 4, quantum paraelectric behaviors are observed. 
The dielectric permittivity follows the so-called Barrett 
formula as a function of temperature, where the diver- 
gence is inhibited by quantum fluctuations and the per- 
mittivity is saturated at low temperatures. Recently, 
photoinduced reflectivity changes are measured^ The 



reflectivity in the energy range where it is sensitive to 
the ionicity change shows a large-amplitude oscillation 
near the quantum phase transition point. 

In this paper, nonequilibrium dynamics are calculated 
after impact force is introduced to a polarization in the 
one-dimensional quantum Blume-Emery-Griffiths model. 
Near the quantum phase transition point, the correlation 
length is long and quantum fluctuations are large because 
the ordering is marginally suppressed. In such a case, 
the impact induces large-amplitude oscillations in polar- 
izations and ionicity. The cooperativity of the induced 
dynamics is measured by calculating the ratios of the 
energies of induced oscillations to the total-energy incre- 
ment. Below, we will show that the dynamics near the 
quantum phase transition exhibit enhanced oscillations 
in amplitudes but their behaviors are uncooperative in 
the present model. 



II. MODEL AND METHOD 

Mixed-stack charge-transfer complexes have different 
types of degrees of freedom, i.e., ionicity and electric po- 
larization, although they are strongly coupled. A global 
change in the ionicity does not alter the symmetry, but 
a finite average of electric polarizations breaks the inver- 
sion symmetry. They can be described by spin-one oper- 
ators, where Sf takes either +1, 0, or — 1. A neutral state 
at site i is represented by Sf=0, while an ionic state with 
positive or negative polarization by Sf = \ or — 1. Here, a 
finite polarization is assumed at each ionic site. Because 
of the spin-Peierls instability, an ionic domain consisting 
of many sites would have dimerization-induced electric 
polarizations. However, a few molecules may not have a 
polarization even if they are ionic. Furthermore, an odd 
number of consecutive ionic (neutral) molecules in the 
neutral (ionic) background cannot be described by such 
a pseudo-spin state. In this sense, the present "sites" 
should be interpreted in a coarse-grained sense. 

Once such coarse graining is accepted, we can em- 
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ploy the Blume-Emery-Griffiths model^ 5 - but one should 
be aware of the following points. If we could integrate 
out electronic degrees of freedom to obtain some effec- 
tive model, it might have some long-range interaction 
that originates from the spin-Peierls instability of one- 
dimensional half- filled correlated electron systems. We 
will not pursue nonequilibrium dynamics specific to par- 
ticular materials, but discuss rather general aspects near 
quantum critical points. Then, we will limit ourselves to 
short-range interactions. We do not distinguish between 
electronic and phonon contributions to electric polariza- 
tions, so that the present "sites" are composites of them. 

The Blume-Emery-Griffiths model is successfully ap- 
plied to TTF-CA. It explains the existence of a triple 
pointy which is experimentally observed*^ The in- 
clusion of interchain Coulomb attraction and inter- 
chain elastic energy qualitatively explains the pressure- 
temperature phase diagram of TTF-CA.— Relaxation 
processes in the neutral-ionic, paraelectric-ferroelectric 
phase transition a 19 ' 20 are discussed on the basis of its 
master equation. 21 

Now, we extend this model to the mixed-stack charge- 
transfer complexes, DMTTF-QBr„Ci4_„, which show 
quantum phase transitions between neutral and ionic 
phases . 12 ' 13 For the spin-one model with short-range in- 
teractions, the way by which quantum natures are in- 
corporated is not unique. Here, we assume that tunnel- 
ing between ionic states with different polarizations is 
achieved always through a neutral state and represented 
by the operator Sf defined below. This quantum version 
is also successful in that the mean-field theory reproduces 
the dielectric permittivity in the neutral phase,— which 
is described by the Barrett formula . 12 ' 13 

This quantum version of the Blume-Emery-Griffiths 
model is defined by 

Hqbeg = SiS- +1 - K y^OSf ) 2 (gf +1 ) 2 

i i 

- P j2(s?) 2 -hj2s? , (i) 

i i 

where Sf = (|1)<1| - | -1><-1 1)< and Sf = (|1)(0| + | 
0)(1 | + | 0)(-l | + | -1)(0 \)i/y/2 with | j)t being the 
eigenstate associated with Sf=j. The periodic bound- 
ary condition is imposed. The parameter J (K) denotes 
the nearest-neighbor dipolar (quadrupolar) interaction, 
P the energy difference between neutral and ionic states, 
and h the neutral-ionic quantum tunneling amplitude. 
Under hydrostatic or chemical pressure, the crystal con- 
tracts and the Madelung energy stabilizes the ionic state 
relative to the neutral state. This would correspond to 
increasing P in the present model. In the P —> oo limit, 
the states 1 0)j are completely suppressed, and the model 
is equivalent to the quantum Ising model (i.e., the spin- 
1/2 transverse Ising model). 

We study dynamics in nonequilibrium conditions by 



introducing the perturbation to Eq. JT]) written as, 

H' = -j2m)S! , (2) 

i 

where Ei(t) = E lmp for i = 1 and < t < t lmp and 
Ei(t) = otherwise. The time-dependent Schrodinger 
equation for the exact many-body wave function is nu- 
merically solved by expanding the exponential evolution 
operator with time slice dt=0.01 to the 15th orde r 10 ' 23 
and by checking the conservation of the norm and of the 
total energy for t > t imp . 

The main purpose here is not necessarily to simulate 
photoinduced dynamics, but to study nonequilibrium dy- 
namics near the quantum critical point from a general 
viewpoint, i.e., coherence and nonlinearity. By the above 
excitation, different polarization states at site i get dif- 
ferent energies. It does not directly transit between dif- 
ferent polarization states. Note that, in itinerant elec- 
tron models, the gauge transformation relates photoex- 
citation through the modulation of site energies via a 
scalar potential to that through the modulation of trans- 
fer integrals via a vector potential^ More importantly, 
Eq. gives a spatially local excitation in contrast to 
realistic photoexcitation, which is regarded as almost a 
uniform excitation. Later, we will discuss how strongly 
coherence is maintained near the quantum phase tran- 
sition point. The impact force introduced above is dis- 
advantageous because the spatial coherence is initially 
deteriorated. Even in such a difficult condition, physi- 
cal quantities will coherently oscillate near the quantum 
phase transition point. 

III. RESULTS 

First of all, the ground-state properties are calculated 
and shown in Fig. [TJ The quantum tunneling between 
any states that are degenerate in the thermodynamic 
limit is inevitable in finite-size calculations and always 
leads to the vanishing expectation value of S z even in 
the case where the symmetry is spontaneously broken 
in this limit. Then, we plot the S z -S z correlation func- 
tion between farthest points, i.e., (S^^ 2 S^) in Fig.QJa). 
It monotonically decreases with increasing h and mono- 
tonically increases with P. Here, the results of differ- 
ent system sizes, N= 8, 10, 12, and 14, are shown for 
each P. The system-size dependence is visible for small 
(Sx/ 2 S]sr)- For large N, it steeply decreases to zero. 
There would be a critical tunneling amplitude h c above 
which (S z N j 2 S^) vanishes in the N — > oo limit. The 
quantity h c increases with P. 

The ionicity ({Sf ) 2 } (at any i) is shown in Fig. QJb). 
For P > —0.4, it monotonically decreases with increas- 
ing h. For any P, it becomes 0.5 in the h — > oo 
limit because the ground state of —h^ ji Sf is (1/2) | 
1)< + (1/V2) | 0) 4 + (1/2) | For P = -0.8, it de- 

creases to a value below 0.5 (0.475) and then gradually in- 
creases to 0.5. With decreasing P, the decrease of ((S z ) 2 ) 
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FIG. 1: (Color online) (a) S z -S z correlation between farthest 
points and (b) expectation value of (S z ) 2 , as a function of 
quantum tunneling amplitude h, with 7=1, K=0.3, and dif- 
ferent P values. For each parameter set, the results of differ- 
ent system sizes, N= 8, 10, 12, and 14, are shown. The curves 
with the steeper slope at the tail correspond to the results of 
the larger size. 



FIG. 2: (Color online) (a) Spatial average of S z and (b) spa- 
tial average of (S z ) 2 , as a function of time t. The quantum 
tunneling amplitude h is varied from 0.6 with a step of 0.1, 
and the corresponding lines are shifted upward by 0.1. Other 
parameters are J=l, A"=0.3, P—~ 1.2, 7V=14, i5i mp =4, and 



=0.3. 



with increasing h becomes steeper. For P = —1.2, it 
abruptly decreases between /i=0.5 and 0.6 and then in- 
creases to 0.5. In the thermodynamic limit, the transition 
from the ferroelectric phase ((Sf) 7= 0) to the paraelec- 
tric phase ((Sf) = 0) would be continuous at least for 
P > —0.8. Meanwhile, the transition seems discontinu- 
ous for P = — 1.2. It indicates the presence of a tricritical 
point between -1.2 < P < -0.8 and 0.5 < h < 1. In 
any case, the dynamical behavior with increasing h above 
the transition point is similar for all P, and the change 
with h is largest for P = —1.2. Then, we will show the 
time-dependent properties for P = —1.2 below. 

In the ground state, the polarization (S 1 ?) is absent 
at any i in the paraelectric phase. We have introduced 
impact force on (5,f) initially and at a particular site 
(i=l), as described by Eq. ([2]). Then, the polarization 
at site 1 grows at first. Owing to the coupling J and 
the periodic boundary condition, the polarizations grow 
successively at sites 2 and N, at sites 3 and N — 1, and 



so on. They oscillate with different phases and interfere 
with each other, so that their spatiotemporal pattern is 
complicated. Because the parameter set used here is in 
the paraelectric phase, their spatial average do not grow 
but oscillate about zero, as shown in Fig.[2ja). Here, the 
lines with different h values are shifted vertically by dif- 
ferent lengths. As the transition point is approached, the 
amplitude becomes large, and the period becomes long 
because the relevant excitation energy decreases. The 
dynamics of polarizations of course affect those of ionic- 
ity. The dynamics of ionicity are more complicated than 
those of polarizations because Eq. @ makes an impact 
directly on the polarization. Despite this, the spatial av- 
erage of ionicity also oscillates, as shown in Fig. ^b), 
though its shape is not as sinusoidal as that of polar- 
izations. There are similarities between the dynamics of 
ionicity and those of polarizations: as the transition point 
is approached, the oscillation amplitude becomes large, 
and the oscillation period becomes long. A detailed com- 
parison will be made later. 
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FIG. 3: (Color online) (a) Spatial average of S z and (b) spa- 
tial average of (S z ) 2 , as a function of time t, with different 
Eimp values. Other parameters are J=l, K—0.3, P— — 1.2, 
h=0.6, N=U, and f imp =0.3. 



In order to see nonlinearity or linearity of the nonequi- 
librium dynamics, its impact-strength dependence is im- 
portant. Figure [3] shows the -Ei mp -dependence of the 
time evolution of the spatially averaged quantities. Here, 
the lines with different Ei mp values are not shifted verti- 
cally. Because the increment in the total energy AE — 
(-Hqbeg) t=t imp - (#qbeg)*=o is approximately propor- 
tional to Ef mp , we vary E? mp with a step of 4. With in- 
creasing -Eimp, the oscillation amplitude becomes large, 
but the oscillation period does not change. The overall 
shape does not vary so much at least in this _Ei, np -range. 
The oscillation amplitude of the spatially averaged ton- 
icity increases almost in proportion to Ef mp [Fig. [3^b)] , 
while that of the spatially averaged polarization increases 
sublinearly with Ef mp [Fig. Efa)] . 

The dependences of these oscillation amplitudes on the 
total-energy increment AE are shown in Fig. 2] for dif- 
ferent h values. Here, the amplitudes are measured at 
the maxima on the first humps (i < 20). The oscillation 
amplitude of the spatially averaged polarization increases 
linearly with the square root of AE [Fig.QJa)], while that 
of the spatially averaged ionicity increases linearly with 



FIG. 4: (Color online) (a) Amplitude of oscillation in spatial 
average of S z , as a function of square root of energy supplied, 
and (b) amplitude of oscillation in spatial average of (S z ) 2 , 
as a function of energy supplied. With decreasing tunnel- 
ing amplitude h toward the quantum phase transition point, 
the slope becomes steep. Other parameters are J=l, K=0.3, 
P=-1.2, iV=14, and t imp =0.3. 



AE [Fig.gjb)]. In Fig.gja), the result for h = 0.6 visibly 
deviates from the linear function, but this is insignificant 
in terms of numerical accuracy. If the amplitude is mea- 
sured at the maximum on the second hump instead, it 
becomes quite a linear function. It should be noted that 
the expectation value of Eq. |T]) contains a term that is 
quadratic with respect to polarizations, and terms that 
are linear and quadratic with respect to ionicity. No term 
is linear with respect to polarizations because of the in- 
version symmetry of the model. Therefore, the result 
shown in Fig.[¥]is a consequence of linear relations among 
the energy supplied by the impulse, the energy flow into 
polarizations, and the energy flow into increased ionic- 
ity. The nonequilibrium dynamics show uncooperative 
behaviors. The impact does not induce a transition into 
the ferroelectric phase. This linearity may be responsible 
for the observed coherence: the phase of the oscillating 
ionicity is maintained even if the impact is made on a 
polarization locally in space and time. If any coopera- 



5 



tive effect were present, the ionicity would not oscillate 
about the initial value and its time evolution would show 
a more complicated behavior. 

What is interesting here is the enhancement of the os- 
cillations near the quantum phase transition point: the 
linear coefficients are increased as the transition point 
is approached. The enhanced oscillation of the spatially 
averaged ionicity is conspicuous. In the quantum para- 
electric phase (more generally, in quantum disordered 
phases), the interaction-induced ordering is suppressed 
by quantum fluctuations. Near the transition point, it 
is marginally suppressed. Then, any small perturbation 
would manifest this fact by showing a large response. 
There are large ferroelectric fluctuations accompanied by 
ionic fluctuations. The correlation length is long near the 
transition point. These facts enhance the coherent os- 
cillations in polarizations and ionicity. These results are 
reminiscent of optical experiments on the perovskite-type 
quantum paraelectric SrTiO a 25 i 26 , which show the pho- 
toinduccd enhancement of the static dielectric permittiv- 
ity near the quantum phase transition point. Although 
SrTiC>3 is a three-dimensional material, significant coop- 
erative behaviors are not observed in the photoinduced 
changes. In the perovskite-type transition-metal oxides 
also, the electron-lattice^ and covalency^ contributions 
on the electronic polarizability are well recognized. 

IV. SUMMARY 

In order to approach the coherence in photoin- 
duced dynamics of quasi-one-dimensional mixed-stack 
charge-transfer complexes near the transition between 
the neutral quantum-paraelectric phase and the ionic 
quantum- ferroelectric phase (within a chain), we the- 
oretically study noncquilibrium dynamics of the one- 
dimensional quantum Blume-Emery-Griffiths model us- 
ing exact many-body wave functions. The time evolu- 
tion of the system is obtained by numerically solving the 



time-dependent Schrodinger equation after impact force 
is introduced on a polarization locally in space and time. 
Even though the impact is made locally, polarizations 
coherently oscillate about zero in the paraelectric phase, 
so that their oscillation survives after taking the spatial 
average. Furthermore, they induce a coherent oscillation 
in the spatially averaged ionicity. 

The induced oscillation is enhanced near the quantum 
phase transition point, so that the oscillation amplitude 
increases with decreasing tunneling amplitude toward the 
quantum phase transition point. The oscillation ampli- 
tude of the spatially averaged polarization increases lin- 
early with the square root of the total-energy increment, 
while that of the spatially averaged ionicity increases lin- 
early with the total-energy increment, implying linear re- 
lations among the energy supplied, the energy flow into 
polarizations, and the energy flow into increased ionic- 
ity. The linear coefficients are large near the quantum 
phase transition point because of large ferroelectric fluc- 
tuations accompanied by ionic fluctuations. Nevertheless 
the linearity implies uncooperative behaviors, so that the 
impact does not induce a transition into the quantum fer- 
roelectric phase. This behavior is reminiscent of photoin- 
duced behaviors of quantum paraelectric SrTiC>3, and it 
is also expected in the photoinduced dynamics in mixed- 
stack charge-transfer complexes near the quantum phase 
transition point. 
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